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ABSTRACT 

We have obtained deep multi-object optical spectra of 49 H II regions in the outer disk of the spiral galaxy 
M83 (=NGC5236) with the FORS2 spectrograph at the Very Large Telescope. The targets span the range 
in galactocentric distance between 0.64 and 2.64 times the R25 isophotal radius (5.4 — 22.3 kpc), and 31 of 
them are located at R >R25, thus belonging to the extreme outer disk of the galaxy, populated by UV com- 
plexes revealed recently by the GALEX satellite. In order to derive the nebular chemical abundances, we apply 
several diagnostics of the oxygen abundance, including R23, [Nll]/[Oll] and the [Olll]A4363 auroral line, 
which was detected in four H II regions. We find that, while inwards of the optical edge the O/H ratio follows 
the radial gradient known from previous investigations, the outer abundance trend flattens out to an approxi- 
mately constant value. The latter varies, according to the adopted diagnostic, between 12 + log(0/H) = 8.2 and 
12 + log(0/H) = 8.6 (i.e. from approximately 1/3 the solar oxygen abundance to nearly the solar value). An 
abrupt discontinuity in the radial oxygen abundance trend is also detected near the optical edge of the disk. 
These results are tentatively linked to the flat gas surface density in the outskirts of the galaxy, the relatively 
unevolved state of the extended disk of M83, and the redistribution of chemically enriched gas following a past 
galaxy encounter. 

Subject headings: galaxies: abundances — galaxies: ISM — galaxies: individual (M83) 



1. INTRODUCTION 

The notion that recent or ongoing star formation activ- 
ity extends well beyond the optical radii of spiral galax- 
ies has developed du ring the last decade, with the detection 
of bo th Hll regions (Ferguso n et alj|l998bt iLelievre & Rov l 
2000) and young B stars dCuillandre et al.1 1200 It iDavidge 
120071) out to approximately twice the 25th magnitude B- 
band isophotal radii (i?2s)- Recent observations with the 
GALEX satellite have also traced the presence of massive 
star formation beyond the R25 radius in the form of UV- 
brigh t clusters and unstructured U V emission (Thil ker et al.l 
120051 120071: iGil de Paz et~aTJl2005l) . Contrary to the paucity 
of galaxies with Hll region detections at large galactocen- 
tric distances, these "extended UV" (XUV) disks are com- 
mon pheno mena, being found in about 1/3 of nearby spi - 
ral galaxies (IZaritsky & Christleinll2007T: IThilker et al.ll2007l) . 
This can be interpreted in terms of the different lifetimes 
of O stars (responsible for the ionization of the Hll re- 
gions, ~ 10 Myr) compared to less massive B stars (pro- 
ducing the bulk of the UV flux, several hundred Myr). Ob- 
servationally, the ages measured for the UV clusters agree 
with this picture, since they range between a few My r and 
- 1 Gyr (IZaritsky & Christleinll2007l: Dong et al.ll2008l) . Al- 
ternative explanations are based on the low star formation 
rate (SFR) and gas density regimes encountered in outer 
spiral disks, which can affect th e upper stellar mass dis- 
tribution in star forming r egions ( Weid ner & Kroupal l200l 
iKrumholz & McKed2008l) . resulting in a lack of massive star 

1 Based on observations collected at the European Southern Observatory, 
Chile, under program 079.B-0438. 

2 Visiting scientist, Institute of Astronomy, University of Cambridge 



production at high galactocentric radii. 

The UV data do not show as pronounced a drop in the SFR 
as the Ha surface brightness profiles do at a radius Ru u (near 
the optical radius), which is usually interpreted as the point 
where the surface density E gas (r) of the gas reaches a crit- 
ical threshold S c (r), below which the gase ous disks of spi- 
ral galaxies become grav itationally stable (iKennicuttl 119891: 
iMartin & Kennicutti2 001). The idea that localized density en- 
hancements [S gas (r) > S c (r)] in subcritical environments can 
produce observable star- forming sites beyond the t hreshold 
radius, as postulated by Martin & Kennicuttj (1200 ll) . is cor- 
roborated by recent modeling that considers either spiral den- 
sity waves prop agating into the extended neutral gas disks 
dBush et al.1 120081) or additional triggering mechanisms for 
star f ormation other than sponta neous gravitational instabil- 
ities felmegree n & Hunterll2006l) . Obser vational support fo r 
this picture has been recently provided bv Dong et al.1 (2008). 
who found that at the locations of the UV-bright clusters in the 
outer disk of M83 the gas density is near the critical value. A 
different interpretation of the discrepant behavior of the Ha 
and UV radial profiles has been pr oposed more recently by 
Pflamm-Altenburg & Kroupa (2008), who invoked the prop- 
erties of clustered star formation and the dependence of the 
maximum stellar mass on the SFR density and cluster mass. 

Hll regions and UV knots located beyond the optical 
edges of spirals n ot only prov i de tra cers of the star for- 
mation processes (IThilker et al.1 120071) and the kinematics 
(IChristlein & Zaritskyll2008h of the outer disk, but also offer 
the opportunity to explore the chemical evolution of a still 
poorly known component of the host galaxies. Infall mod- 
els of galaxy formation predict t hat spiral disks build up via 
accretion by growing inside-out (Matteucci & Francois 1989; 
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iBoissier &^ Prantzos 1999), so that the outermost parts of the 
disk should possess low metallicities and could possibly be 
still subject to active phases of infall. Consequently, ob- 
servational data on galactic abundance gradients along ex- 
tended radial baselines can provide s trong constraints for 
chemical evolution models of galaxies (Chiappini et al JI200 lb 
ICescutti et af]|2007h . 

The few studies of the chemical composition of H II regions 
at large galactocentric distances that are available in the liter- 
ature suggest that the extended disks are relatively unevolved 
systems. In their spectroscopic work on th e late-type spi- 
rals NG C 628, NGC 1058 and NGC 6946, IFerguson et aH 
(1998a) found that the abundance of oxygen (a proxy for 
the overall metallicity in ionized nebulae) decreases in a 
log-linear fashion from 12 + log(0/H) = 8.9 near the center 
to 12 + log(O/H) = 8.0 for the outermost (1.5-2. 0R 25 ) targets 
[i.e. from nearly twice the solar abund ance down to 1/5 sol ar, 
if 12 + log(O/H) =8.66, the value bv lAsplund eTal]|2005L is 
adopted]. The N/O ratio they measured in the outermost re- 
gions is approximately 20-25% of the solar value. The major 
limitation of this early study is represented by the fact that in 
the three target galaxies only a total of nine H II regions lo- 
cated at R > R25 were analyzed, so that general conclusions 
about the shape of the radial abundance gradients and the 
scatter in metallicity become very uncertain. The low oxy- 
gen ab undances measured in outer disks by IFerguson et al.1 
(11998a!) agree, however, with the values o btained near the op- 
tical edges of a few spiral s (e.g. M101: iGarnett & Kennicuttl 
11994b Ivan Zee et al.lll998b iKennicutt et al.ll2003l) . More re- 
cently, Gil de P az et al.l (120071 = G07) have reported similarly 
low metallicities (10-20% of the solar value for oxygen) for 
H II regions in the extended disks of M83 and NGC 4625, but 
with a high N/O ratio, close to solar, in the case of M83. 

Measuring chemical abundances in extended spiral disks 
is technically challenging. The ty pical Ha luminosity of 
the outermost targets inve stigated by IFerguson et al.1 {1998a) 
and iLelievre & Royl (l2000h is L(Ha) < 10 38 ergs _1 , which 
corresponds to the ionizing output of a small number of O 
stars. These luminosities translate into extremely faint fluxes 
in the emission lines that are used as metallicity diagnos- 
tics. In particular, the detection of the auroral lines (such 
as [O III] A4363) used in the classical chemical abundance 
analys is of ionized nebulae to deri ve the electron temperature 
(e.g. lPagellll997b ISta siriska 2007) would become extremely 
difficult in all but the brightest outer disk H II regions in galax- 
ies at distances of just a few Mpc. Statistical indicators of 
gas abundances, based on stronger nebular lines, can be em- 
ployed as an alternative. For instance, taking advantage of 
the presence of both singly and doubly ionized oxygen line 
transit ions in the optical spectra of H II regions, Pagel et al. 
(1 19791) introduced the metallicity i ndicator R^ = (\Q II] A3 727 
+ [O in] AA4959, 5007)/H/3 (both IFerguson et alj|l998al and 
G07 base their outer disk chemical abundances on R23). As 
already explained in detail in a number of papers, the calibra- 
tion of these indexes in terms of oxygen abundance is prob- 
lematic, and c an lead to still poorly understood systematic 
discrepancies dPerez-Montero & Dfazl 120051 iBresolinl 12008b 
iKewlev & E llison 20081). 

In this work we look at the chemical abundance proper- 
ties, in particular the radial oxygen abundance gradient and 
the N/O abundance ratio, of the outer disk of M 83, the pro- 
totypical XUV disk galaxy dThilker et al.1 120051) . While the 
number of UV complexes with associated H II emission is 
only a fraction of the total of more than 100, there are still 



several dozen faint Hll region candidates beyond the opti- 
cal edge of the galaxy that are visible in deeply exposed Ha 
frames. Within the context of a project that aims at obtain- 
ing chemical abundance properties in extended spiral disks, 
we have targeted these emission-line objects with a multi- 
slit spectroscopic program. While a similar project was pub- 
lished by G07 at the time we started to acquire our data, we 
overcome some of the deficiencies encountered in their work, 
by enlarging the sample of outer disk H II regions (> 30 tar- 
gets at R > R25), and by obtaining deeper spectra. In order 
to facilitate the comparison with inner disk nebular abun- 
dance s in M83 dBresolin & Kennicuttl l2002t iBresolin etafl 
we have also included in our sample a number of H II 
regions that are located inside the Run radius (taken from 
iMartin & Kennicuttl 12001 1: 5.1 = 6 .7 kpc at the adopted dis- 
tance of 4.5 Mpc. lThim et al.ll2003l) . Our principal goal at the 
onset of the project was to verify whether in the prototypical 
XUV disk galaxy the oxygen abundance decreases steadily 
with galactocentric distance, down to values comparable to 
those found in previous investigations (10-20% of the solar 
value). 

Here we present the results of our study of the outer disk 
of M83. The observational material and the analysis of the 
multi-object spectroscopy are presented in §2. In §3 we 
carry out the chemical composition analysis, using a vari- 
ety of metallicity indicators. Our discussion concerning the 
observed oxygen abundance gradient in M83, from the inner 
part of the disk to the outermost H II region (at 2.6 ^25), and 
the interpretation of the results, are presented in §4. Lastly, 
we summarize our findings in §5. 

2. OBSERVATIONS AND DATA REDUCTION 

Candidate H II regions in the outer disk of M83 were se- 
lected from deep Ha images obtained at the 90-inch Bok tele- 
scope of the Steward Observatory with the 90Prime camera. 
Several of these Hll regions have a counterpar t in the list 
of UV -emitting clusters detected by GALEX (Thilk er" ;t alJ 
12005b G07). The azimuthal distribution of the candidates 
is not uniform, since, as shown in the GALEX images, the 
brightest star-forming sites in the outer disk of M83 are orga- 
nized in discrete filamentary structures that extend outwards 
from the edge of the optical disk of the galaxy, in general cor- 
responding with enhancements in the radio H I emission. 

Optical spectra were acquired with FORS2 at the Very 
Large Telescope of the European Southern Observatory on 
Paranal. Objects distributed in four non-overlapping 6.'8 x 
6. '8 fields were observed in service mode between May 10, 
2007 and June 19, 2007. One of these fields covers part of the 
northern extended disk of the galaxy, where a faint UV fila- 
ment, resembling an external spiral arm, extends to the north 
of the main disk of the galaxy. Two fields are located to the 
south (one of them partially overlapping with the bright op- 
tical disk), while the last field is located to the south-west of 
the optical disk (Fig. [TJ. The slits of the multi-object masks 
used for the observations included, in addition to some promi- 
nent UV-emitting knots, a number of additional emission-line 
objects identified on Ha pre-imaging frames. In filling the 
slit masks priority was given to the outermost H II region can- 
didates. The spectra were obtained through 1.0 arcsec-wide 
slitlets with the 600B grism ('blue' setting, central wavelength 
4650 A, 6 x 970 s exposures) and with the 1200R grism ('red' 
setting, central wavelength 6500 A, 3 x 335 s exposures), pro- 
viding an almost full spectral coverage from approximately 
3500 A to 7200 A for most of the targets (the exact spectral 
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TABLE 1 
HII REGION SAMPLE 



ID 


R.A. 


Decl. 


R/R 25 


Notes 




(J2000.0) 


(J2000.0) 






(1) 


(2) 


(3) 


(4) 


(5) 




FIG. 1. — Distribution of the four fields in which the multi-object spec- 
troscopy has been carried out. Each field is 6. 8 X 6' 8 wide. The image is a 
composite of GALEX and Ha data, with UV-bright knots appearing in blue, 
and H II regions in red. 



coverage depending on the location of the objects within the 
field). With these setups, the FWHM spectral resolution is 
4.5 A and 2.4 A for the blue and red spectra, respectively. The 
red and blue observations were carried out in consecutive ex- 
posures, thus minimizing possible differences in target cov- 
erage arising from variations in the seeing. Spectra for the 
southern field were not obtained in the same night, but the im- 
age quality during the observations was similar. Typical see- 
ing conditions during the observations were ~ 0.8 arcsec. The 
use of the FORS atmospheric dispersion corrector dAvila et all 
1 19971) allowed us to avoid the adverse effects of differential re- 
fraction on the emission line ratios and the derived chemical 
abundances. Observations of the standard star LTT 4816 have 
been obtained to flux-calibrate the H II region spectra. 

For the data reduction we relied on the flat-fielded and 
wavelength calibrated products of the EsoRex pipeline (ver- 
sion 3.1.1). We then obtained flux-calibrated spectra, clean 
of cosmic rays, using common IRAF 3 routines. Since the 
standard star spectra were obtained through a single slit near 
the center of the detector, we investigated possible variations 
in the response function across the CCD by comparing the 
wavelength dependence of flat fields obtained from slitlets 
distributed across the FORS2 field. Once scaled to a com- 
mon mean flux, the resulting rms dispersion of the flat field 
spectra is ~ 0.03 mag, which is also our adopted uncertainty 
for the flux calibration. 

Emission-line fluxes were measured for 49 H II regions be- 

3 IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in As- 
tronomy, Inc., under cooperative agreement with the National Science Foun- 
dation. 



XUV01 



1 13 35 56.71 -29 56 37.3 2.33 

2 13 36 12.45 -29 57 06.6 1.84 

3 13 36 13.27 -29 52 44.6 1.68 

4 13 36 19.25 -29 53 43.8 1.47 

5 13 36 28.58 -29 56 06.0 1.28 

6 13 36 29.23 -29 54 12.2 1.15 

7 13 36 29.30 -29 54 34.1 1.16 

8 13 36 31.08 -29 55 42.0 1.17 

9 13 36 32.53 -29 54 13.0 1.04 

10 13 36 46.55 -29 46 40.5 1.03 

11 13 36 46.75 -29 44 42.0 1.32 

12 13 36 47.71 -29 46 28.1 1.04 

13 13 36 49.08 -29 56 02.9 0.76 

14 13 36 50.07 -29 55 52.7 0.72 

15 13 36 50.25 -29 42 28.5 1.62 

16 13 36 51.10 -29 56 13.1 0.75 

17 13 36 53.09 -29 55 50.0 0.67 

18 13 36 54.24 -29 47 13.8 0.82 

19 13 36 54.79 -29 47 07.2 0.83 

20 13 36 55.17 -30 05 47.4 2.25 

21 13 36 55.35 -29 47 32.1 0.76 

22 13 36 55.64 -29 47 37.6 0.74 

23 13 36 55.69 -29 55 46.4 0.64 

24 13 36 56.79 -29 55 52.2 0.65 

25 13 36 57.23 -30 08 12.0 2.64 

26 13 36 57.59 -29 41 07.7 1.77 

27 13 36 58.07 -30 07 40.2 2.56 

28 13 36 58.17 -30 07 17.1 2.49 

29 13 36 58.52 -29 59 24.5 1.21 

30 13 36 58.64 -30 06 00.5 2.29 

31 13 36 59.35 -29 41 18.7 1.73 

32 13 36 59.62 -30 06 21.0 2.34 

33 13 37 00.29 -29 41 24.0 1.72 

34 13 37 00.39 -29 47 45.9 0.68 

35 13 37 00.49 -29 57 23.2 0.89 

36 13 37 01.19 -29 57 10.8 0.86 

37 13 37 01.69 -30 00 54.3 1.55 

38 13 37 04.47 -29 46 55.9 0.81 

39 13 37 04.99 -29 59 45.9 1.30 

40 13 37 05.82 -30 00 11.2 1.37 

41 13 37 06.76 -29 59 57.0 1.34 

42 13 37 07.13 -29 56 47.6 0.84 

43 13 37 08.20 -29 59 19.6 1.25 

44 13 37 08.93 -29 56 39.8 0.84 

45 13 37 09.74 -29 59 04.5 1.23 

46 13 37 10.76 -29 55 45.4 0.74 

47 13 37 12.36 -29 57 52.9 1.08 

48 13 37 16.42 -30 00 47.6 1.59 

49 13 37 19.35 -29 57 40.3 1.19 



XUV 13 
XUV02 



XUV 06 

XUV 09 

XUV 08 
XUV 22 
XUV 11 



NOTE. — Units of right ascension are horns, minutes and seconds, and units of 
declination are degrees, arcminutes and arcseconds. Col. (1): Hit region identification. 
Col. (2): Right Ascension. Col. (3): Declination. Col. (4): Galactocentric distance in 
units of R25 = 6.44 arcmin. Col. (5): Identification from G07. 



longing to M83. Their coordinates are presented in Table [T] 
together with the galact ocentric distances R/R25 (adopting 
R25 = 6.44 arcmin from Ide Vaucouleurs et"ai1ll99lL an incli- 
nation angle i = 24° from lHeidmann et alJI 19721 and a value 
for th e position angle of the major axis 9 = 45° from lComtel 
I198U we did not accoun t for the warp of the disk at large 
radii. iRogstad et all 1 19741) . and the identification of the ob- 
jects in common with the study by G07. As Table Q] shows, 
our sample covers H II regions from the outer third of the main 
disk of star formation of M83, out to 2.6 times the isophotal 
radius (/?2s)- A small number of targets, although detected in 
the Ha line, were too faint to allow a measurement of addi- 
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tional lines in their spectra, and were therefore discarded from 
the subsequent analysis. We also removed from the final sam- 
ple a number of background line-emitting galaxies (listed in 
the Appendix). 

For most of our targets we were able to measure the in- 
tensities of the strongest recombination lines of the hydro- 
gen Balmer series, and the collisionally excited metal lines 
[On]A3727, [O m] AA4959, 5007, [Nn]AA6548, 6583, and 
[S n] AA6717, 6731. In a handful of cases the signal-to-noise 
ratio of the spectra was sufficient to detect some of the fainter 
metal lines, such as [Ne III] A3869 and [Arm] A7135. In four 
H II regions we also measured the [O III] A4363 auroral line, 
which, in combination with [O III] AA4959, 5007, provides a 
direct determination of the electron temperature of the ionized 
gas (as described in §|4). Due to the moderate redshift of the 
galaxy (v = 513 kms" 1 , Koribalski et al. 2004), the measure- 
ment of the strongest helium recombination line, He I A5876, 
was hampered by the contamination of the sodium telluric D2 
line at 5890 A. The other helium lines were generally too faint 
to be measurable. 

As a rule, the blue and red spectra have no emission lines 
in common. We then decided to refer the strengths of the 
emission lines in the blue wavelength range (in particular the 
[O II] and [O III] lines) to the strength of H/3, while in the 
red we express the line strengths of [N II] and [S n] relative 
to Ha, and assume the case B ratio Ha/H/3 = 2.86 (valid at 
T e = 10 4 K) in order to match the two sets of lines (this renor- 
malization was done after the correction for reddening was 
applied to the individual line fluxes). Since, with the excep- 
tion of [O II] A3727, all the key forbidden lines used for the 
metallicity determination are close to their reference Balmer 
lines, the diagnostic ratios are fairly insensitive to reddening. 
The interstellar reddening was determined from the observed 
H7/H/3 ratio (these two Balmer lines have been detected in 
all 49 H II regions) and, when allowed by the signal-to-noise 
ratio, the HSfHf3 ratio, compared to the theoretical values for 
case B at 10 4 K. The resulting reddening-corrected line fluxes, 
expressed in units of H/3 = 100, are presented in Table|2] For 
the [O III] and [N II] doublets we provide the fluxes of the 
strongest line of the pair ( A5007 and A6583, respectively). 
For [S II] we provide the sum of the A6717 and A6731 fluxes. 
The data for [On]A3727 and [Arm] A7 135 are sometimes 
missing due to incomplete wavelength coverage. The errors 
quoted in Table |2]include, besides the statistical errors and an 
estimate of the flux-calibration errors, also the uncertainty in 
the reddening correction, which could be significant in partic- 
ular for the [O II] A3727 line. In general, the extinction values 
(shown in column 8) are small, with an average c(H/3) = 0.08 
for/?25> 1, compared to an average c(H/3) = 0.17 for 7?25 < 1. 

The H/3 line flux in column (7) is the flux measured within 
the 1-arcsec slits, corrected for extinction, and likely underes- 
timates the total nebular flux, due to slit slosses. As already 
pointed out by G07, these faint nebular fluxes correspond to 
the ionizing output of one late-O star or, at most, a few mas- 
sive stars. Fig. [2] compares the histogram of the extinction- 
corrected number of ionizing photons Qo corresponding to the 
H/3 fluxes in Table [2] (for a distance of 4.5 Mpc) to that ob- 
tained from two samp l es of i nner disk H II regions stud ied by 
iBresolin & Kennicuttl d2002l) and iBresolin et all d2005l) . As 
the histograms show, the median value of our outer disk sam- 
ple is log Qo = 48. 17 (approxim ately equivalent to th e ionizing 
output of a single 08.5 V star, iMartins et"ai1l2005l) . while in 
the giant H II regions of the inner disk the median value is 
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FIG. 2. — Histograms showing the distribution of the extinction-corrected 
number of ionizing photons Qo (in s ) for the Hll regions in the outer disk 
of M83 (green shaded area, left) and for the samples of inn er disk regions 
from Bresolin & Kennicutt 12002) and Bresolin et al. 12005, yellow shaded 
area, right). The median values of the two distributions are indicated by the 
red marks. 



two orders of magnitude larger, log Qo = 50.21 (~130 08.5 V 
stars, or four 03 V stars). It should be noted that the aper- 
ture corrections are likely to be larger for the supergiant H II 
regions that make up the inner disk sample compared to the 
more compact outer disk regions, which would increase the 
offset seen in Fig.|2]between the two samples. 

A quick characterization of the H II regions in the outer disk 
of M83 can be obtained by plotting the traditional [N Il]/Ha 
vs. [Oi hI/H/3 and [S ill/Ha vs. [Oml/H/3 diagnostic di - 
agrams (iBaldwin et all I198U IVelileux & Osterbrockl 119871) . 
Fig.[3]shows that the excitation sequence formed by the outer 
disk Hll regions merges into the sequence outlined by the 
inner disk H II regions, which are found at lower excitation 
(smaller [O m]/H/3 ratio), likely as a result of the higher metal 
content near the center of the galaxy. Two objects, number 5 
and 36 in our list, deviate from the main sequences of points 
in Fig. [3] and lie above the curves that define the upper bound- 
aries of the loci occ upied by star-forming regions (taken from 
iKewley et al.l2006h . The ionization structure of these two H II 
regions appears to differ from that of the remaining objects, 
and in particular the high [N ll]/Ho; and [S ll]/Ho! ratios mea- 
sured in object 36 resemble the signatures of shock heating 
typically observed in supernova remnants. This object was 
in fact already identified as a supernova rem nant candidate 
on the basis of the strong [S Il]/Ha ratio by iBlair & Lond 
(2004, their object BL34). Its emission line ratios are con- 
sistent with this conclusion if one compares the location of 
supernova remnants and photoionized ne bulae in a [S ll]/Ho; 
us. [Nn]/Ha plot (ISabbadin et al.lll977l) . We will remove 
these two targets (5 and 36) from the remainder of the discus- 
sion and from all of the following diagrams. 

We have compared our fluxes in Table [2] with the values 
published by G07 for the eight targets in common. The agree- 
ment is satisfactory for the [Nil], [S II] and [Om] fluxes in 
most cases. An obvious discrepancy is represented by our 
target #30 (XUV 02 in G07), with our [O m] A5007 flux be- 
ing more than 6 x smaller. Although the coordinates of our 
target match, within the uncertainties, those by G07, we be- 
lieve that the slits did not actually include the same object 
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TABLE 2 
Reddening-corrected fluxes 



ID 


[Oil] 


[Om] 


[Nil] 


[Sir] 


[Arm] 


F(H/3) 


c(H/3) 




3727 


5007 


6583 


6717+6731 


7135 


(erg s~' crrT 2 ) 


(mag) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 



1 370 ± 20 146 ± 7 36 ± 3 

2 472 ± 59 197 ± 13 63 ± 12 

3 ■ • • 29 ± 4 60 ± 12 

4 ■•• 162 ±7 44 ±3 

5 716 ± 86 593 ±58 94 ± 18 

6 ■•• 21 ±2 96 ±8 

7 ■ ■ • 99 ± 5 70 ± 7 

8 265 ± 26 248 ± 12 66 ± 7 

9 ■■• 11 ±1 90 ±7 

10 338 ±21 121 ± 6 50 ±5 

11 509 ± 78 36 ±5 51 ±32 

12 639 ± 109 57 ± 7 88 ± 15 

13 214 ± 19 108 ± 6 107 ± 10 

14 113 ±7 14 ±1 96 ±5 

15 440 ± 59 73 ± 5 58 ± 15 

16 94 ± 6 I) : I 85 ± 6 

17 211 ± 12 44 ±2 138 ±9 

18 206 ± 25 22 ±4 114 ±14 

19 216 ± 22 58 ±4 110 ±9 

20 348 ± 18 287 ± 13 37 ± 3 

21 184 ± 17 42 ±3 118 ±12 

22 185 ±20 17 ±3 126 ± 10 

23 128 ±7 50 ±2 92 ± 6 

24 76 ± 7 ± 1 79 ± 6 

25 485 ± 38 123 ± 6 49 ± 8 

26 464 ± 35 67 ± 4 46 ± 5 

27 461 ± 29 302 ± 14 47 ± 5 

28 381 ±77 0±6 49 ± 25 

29 254 ± 14 192 ± 9 59 ± 3 

30 382 ±22 110 ± 5 55 ± 5 

31 389 ±21 105 ± 5 46 ± 3 

32 716 ± 104 0±4 57 ± 17 

33 342 ±29 112 ± 6 54 ± 12 

34 213 ±14 19 ±1 121 ±7 

35 101 ± 8 ± 1 99 ± 7 

36 1071 ± 82 269 ± 13 324 ± 23 

37 ■■■ 23 ± 3 85 ±18 

38 88 ±21 0±3 108 ±13 

39 257 ± 13 258 ± 12 37 ± 2 

40 352 ±19 192 ±9 43 ± 3 

41 159 ± 29 340 ±26 53 ±11 

42 198 ± 10 60 ± 3 121 ± 5 

43 336 ± 18 274 ± 12 39 ± 3 

44 193 ±15 13 ±2 127 ± 12 

45 539 ± 37 25 ± 2 59 ± 3 

46 97 ±5 6 : I 111 ±6 

47 280 ± 36 43 ±4 114 ±14 

48 ■■■ 121 ± 6 62 ±6 

49 399 ± 29 33 ± 2 105 ± 8 



28 ±2 




1.7 xlO- 15 


0.00 


78 ± 15 


0.0 ± 15.8 


3.0 XlO- 16 


0.08 


69 ± 13 


0.0 ± 11.3 


1.1 xlO- 16 


0.17 


53 ±3 


9.6 ±2.1 


7.7 XlO- 16 


0.05 


96 ± 18 


52.8 ± 16.8 


5.0 XlO" 17 


0.00 


76 ±6 


0.0 ± 3.2 


2.4 XlO" 16 


0.00 


57 ±5 


12.2 ± 5.1 


5.1 XlO" 16 


0.11 


57 ±7 


0.0 ± 4.4 


3.7 XlO" 16 


0.32 


80 ±7 


0.0 ± 4.9 


2.6 XlO" 16 


0.00 


23 ±6 




4.4 XlO" 16 


0.00 


81 ±29 




6.5 XlO- 17 


0.00 


0± 13 




2.8 XlO" 16 


0.39 


57 ±7 




1.8 XlO" 16 


0.20 


29 ±3 




1.4 xlO- 15 


0.14 


30 ± 12 




7.3 XlO" 17 


0.00 


42 ±5 




2.4 XlO- 16 


0.00 


84 ±6 




2.9 xlO" 16 


0.00 


60 ±7 




1.5 xlO- 16 


0.15 


44 ±4 


0.0 ± 2.7 


4.0 X 10- 16 


0.35 


30 ±2 


10.8 ± 1.9 


1.6 xlO" 15 


0.00 


39 ±5 


0.0 ±3.8 


2.0 xlO- 16 


0.00 


83 ±6 


0.0 ± 2.8 


7.9 XlO" 16 


0.47 


40 ±3 




6.2 xlO" 16 


0.04 


50 ±5 




4.0 x 10- 16 


0.12 


91 ± 11 


0.0 ±7.8 


4.6 xlO" 16 


0.11 


43 ±5 


0.0 ± 2.9 


3.1 xlO" 16 


0.00 


65 ±6 


0.0 ± 4.7 


6.1 xlO" 16 


0.05 


0±37 


0.0 ± 27.8 


5.1 xlO" 17 


0.00 


30 ± 1 


9.6 ±0.6 


5.6 xlO- 15 


0.00 


38 ±5 


0.0 ± 4.7 


7.5 xlO -16 


0.12 


44 ± 2 


6.6 ± 1.2 


7.2 xlO" 16 


0.00 


111 ± 27 


0.0 ± 2.0 


1.0 xlO- 16 


0.00 


103 ± 11 


0.0 ± 5.7 


1.9 xlO- 16 


0.00 


62 ±4 


0.0 ± 1.4 


1.3 xlO" 15 


0.23 


41 ±4 




2.8 xlO- 16 


0.05 


224 ± 12 




3.0 xlO- 16 


0.33 


110 ±22 


0.0 ± 12.9 


8.7 XlO" 17 


0.05 


82 ±9 


0.0 ± 5.2 


2.2 XlO" 16 


0.42 


32 ± 1 


9.4 ± 0.7 


1.3 XlO" 15 


0.00 


38 ±2 


12.9 ± 1.8 


3.8 XlO- 16 


0.09 


65 ± 14 


0.0 ± 9.8 


7.3 XlO" 17 


0.00 


36 ± 1 




2.4 XlO" 15 


0.31 


30 ±2 


8.5 ± 1.6 


5.4 XlO" 16 


0.00 


58 ±7 




3.1 XlO" 16 


0.00 


65 ±4 


0.0 ± 2.0 


5.5 XlO" 16 


0.60 


51 ±2 




1.3 XlO" 15 


0.05 


70 ± 11 


0.0 ± 8.7 


8.7 XlO- 17 


0.09 


51 ±6 


0.0 ±5.3 


2.7 XlO" 16 


0.00 


64 ±6 




3.0 XlO- 16 


0.19 



NOTE. — Line fluxes in columns (2) and (3) are in units of H/3 = 100, those in columns (4)-(6) are in units of Ha = 286. The extinction correction has been applied to F(H/3) and the 
error on c(H/3) has been propagated to the line fluxes in columns (2)-(6). 



or the same portion of it. In all cases we measure stronger 
[O II] A3727 fluxes (5 of the 8 objects in common have a 3727 
flux in G07), despite the fact that we derive smaller extinction 
values (0.06 vs. 0.16 for the average c(H/3) of the objects in 
common). The worst case (factor of 4 discrepancy) is rep- 
resented by #42 (XUV 22 in G07). For this target the small 
signal-to-noise ratio in the G07 spectrum is likely the reason. 
It is well known that [O II] A3727 line fluxes are the most un- 
certain in nebular spectra, being the most affected by errors 
in the flux calibration and in the interstellar reddening correc- 



tion. While it is difficult to assess the accuracy of the line 
ratios, we have verified in the few cases were this was possi- 
ble that the strengths of high order Balmer lines (H9, H10 and 
HI 1) relative to H/3 agree with the theoretical Balmer decre- 
ment. This provides some confidence about the reliability of 
our fluxes for the nearby [O II] A3727 line. 

3. NEBULAR ABUNDANCES 
3.1. Electron temperatures and direct abundances 
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FIG. 3. — Emission-line diagnostic diagrams for the M 83 outer disk sample 
(re d dots) and for the inn er disk H II regions of Bresolin & Kennicuti] J2002I) 
and Bresolin et al. 12005, open light blue circles). Th e curves repre s ent th e 
upper boundaries for photoionized nebulae defined by Kewley et al] 120061) . 
The two outlying regions, #5 and #36, are identified. 

In order to derive reliable chemical abundances of ion- 
ized nebulae it is necessary to have a good knowledge of the 
physical conditions of the gas, in particular of the electron 
temperature T e , because the line emissivities depend strongly 
on it. Nebular electron temperatures can be obtained from 
lines of the same ions that originate from different excita- 
tion levels, such as the auroral [O III] A4363 and the nebular 
[Om] AA4959, 5007 lines. In the case of high-excitation (low 
metallicity) extragalactic H II regions and planetary nebulae 
the [O III] A4363 line is often detected, but it becomes unob- 
servable as the cooling of the gas becomes efficient at high 
metallicity, or whenever the objects are faint. 

Because of the intrinsically weak levels of line emission 
from the H II region population in the outer disk of M83, in 
which the ionizing sources appear to be single stars or very 
small clusters, it is not surprising that the A4363 line remains 
undetected for the vast majority of the objects. The chemical 
abundance analysis must then be based on statistical meth- 
ods, without knowledge of the electron temperature (as dis- 
cussed in the next section). However, we obtained significant 
[O III] A4363 detections in four H II regions, which allowed 
us to derive their oxygen abundance via the direct method 
based on the knowledge of T e . In Table [3] we present their 
identification and the measured A4363 fluxes. Using the rou- 
tines in the nebular package of IRAF we have derived the 
electron temperature (column 3), the total oxygen abundance 
12 + log(0/H) (column 4), and the N/O abundance ratio (col- 
umn 5), assuming that N/O = N + /0 + . 

While we postpone a full discussion of the chemical com- 
positions to the next section, we briefly note that: (a) one of 
the outermost H II regions in our sample, #20 at R/R25 = 2.25, 
has an abundance 12 + log(0/H) = 8.17, equivalent to 1/3 of 
the solar value. Two of the remaining objects, #39 and 
#43, located between 1.25 and 1.30 isophotal radii, have 
similar oxygen abundances, while the remaining #29 ap- 
pears to be even more metal-rich (half-solar abundance), 
despite being located well beyond the isophotal radius at 
R/R25 =1.21. The results from this limited number of A4363 
detections suggest that the outer disk of M83 is, perhaps 



TABLE 3 
Direct method abundances 



ID 


[Om] 




12 + log(0/H) 


log(N/0) 




4363 


(K) 






(1) 


(2) 


(3) 


(4) 


(5) 


20 . 


2.8 ± 0.3 


11550 ±660 


8.17 ±0.07 


-1.18 ±0.08 


29 . 


0.7 ± 0.1 


8730 ± 470 


8.41 ± 0.08 


-0.98 ± 0.09 


39 . 


2.7 ± 0.3 


11810 ± 640 


8.05 ± 0.07 


-1.03 ±0.08 


43 . 


2.9 ± 0.6 


11895 ± 1080 


8.12 ±0.12 


-1.14 ±0.14 



NOTE. — The reddening-corrected intensity of [O III] A4363 is in units of H/3 = 100. 



somewhat surprisingly, not as metal-poor as expected. G07 
assigned an abundance 12 + log(0/H) = 7.86 to region #20 
(their NGC 5236:XUV01) based on the value of the abun- 
dance indicator R23. Our metallicity is a factor of two larger. 
(b) the N/O ratio ranges approximately b etween 50% and 80% 
of the solar value riog(N/O) =-0.88, lAsplund et all 12005m . 
As we will show in the next section, these N/O ratios are sim- 
ilar to those found in extragalactic H II regions of comparable 
O/H abundances. 

3.2. Strong-line abundances 

In the absence of direct indicators of the electron tempera- 
ture such as the auroral lines, the chemical analysis must be 
carried out via methods that statistically determine an oxy- 
gen abundance from the ratios of strong nebular lines. These 
strong-line indicators and their calibration in terms of oxy- 
gen abundance have been discussed at length in recent years, 
due especially to the necessity of resorting to these methods 
for the measurement of abundances in star-forming galaxies 
at intermediate and high redshift. In order to obtain a pic- 
ture of the chemical composition in the outer disk of M83, 
we will consider a number of strong-line indicators. It is 
well-known that different indicators provide different solu- 
tions to the chemical composition of star-forming galaxies 
and Hll regions. Some of the methods only use oxygen 
lines to estimate oxygen abundances, while others use lines of 
other elements (e.g. nitrogen) as proxies of the oxygen abun- 
dance. Moreover, the same strong-line indicator can be cali- 
brated empirically from nebulae where a T^-based abundance 
is available, or from grids of photoionization models. Once 
again, systematic differences are found, with models provid- 
ing higher abundances (u p to 0.7 dex) than r e -based methods 
(Kewley &Ellison||2008). Some indicators are strongly de- 
pendent on the ionization parameter, while others are almost 
insensitive to it. By considering different strong-line indica- 
tors we attempt to provide robust results concerning the trends 
of chemical composition in the extended disk of M83, with 
the warning that a significantly more detailed picture is be- 
yond our current capabilities. One potential difficulty in the 
use of strong-line indicators for our sample of low-luminosity 
Hll regions is that the available calibrations are based on 
more luminous giant and supergiant H II regions. However, 
tests of strong-line diagnostics in Galactic nebulae ionized by 
one or a few stars show no evidence for systematic effects 
dKennicutt et al.l2000HOe~v & Shieldsl2000T) . As a further em- 
pirical test, we have applied one of the diagnostics described 
below ([N ll]/[0 II]) to a number of well-studie d H II regions 
in the Milky Way, including the Orion nebul a (Esteba n et alj 
120041) . M8, M17 and others (see lGarcfa-Roias et al.ll2007l and 
references therein). Admittedly, for the sample considered 
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12 + log(0/H) only ranges between 8.39 and 8.56. With the 
particular [N ll]/[0 II] calibration adopted (lBresohrj2007l ob- 
tained from extragalactic giant H II regions) we were able to 
recover the published T f -based O/H abundances within 0.12 
dex in all cases, with a mean difference of 0.04 dex and a stan- 
dard deviation of 0.05 dex. The smallness of the latter value 
indicates that the adopted strong-line diagnostic is quite reli- 
able also when applied to lower luminosity nebulae, at least 
near the O/H range of the Galactic sample. 

The strong-line indicators considered in this work are: 

(a) R 2 3=([Oll]A37 27 + [O m] AA4959,5007)/H /3. Despite 
its drawbacks (e.g. Perez-Monte ro & Dfazl 12005), this clas- 
sic indicator is widely used. One of the main difficul- 
ties in its adoption is related to the necessity of locat- 
ing which of two branches (upper and lower) a given 
Hll region belongs to, since R23 is degenerate. This 
is usually accomplished by simultaneously looking at 
line ratios that increase monotonically with oxygen abun- 
dance, such as [N Il]/H a or [Nll]/[Oll]. According to 
iKewlev & Ellison! (|2008), an object belongs to the upper 
branch of R 23 if log([N 11] A6583/[0 11] A3727) > -1 .2 and/or 
log([N II] A6583/Ha) > -1 . 1 [the lowest values in our sample 
are log([Nll]/[Oll]) = -l.l, log([Nll]/Ha) = -0.90]. Accord- 
ing to this scheme, all of the H II regions in our sample belong 
to the R23 upper branch [12 + log(0/H) = > 8.4 in the abun- 
dance scale calibrated via photoionization models]. For the 
derivation of oxygen abundanc es we adopt the upper branch 
ca libration of McGaugh (1991), in the analytical form given 
bv lKobulnickv et al.1 dl999T ~ 

(b) N2 = [N III A6583 /Ha, as empirically calibrated by 
iPettini &P agel (2004). This indicator, like the remaining 
ones considered by us, is a monotonic function of oxygen 
abundance, thus avoiding the degeneracy problem of R23. 

(c) 03N2 = logr(rOin1 A5007/H/3)/ (rNii] A6583/Hq)], also 
calibrated by Petti ni & Pagel (120041) . 

(d) r Nii1A 6 583/rO il1A3727, which, as shown by 
Dopit a eUri] (|2000), is highly insensitive to the ioniza- 
tion parameter. W e adopt here two d i fferent calibrations, the 
theoretical one by IKewle v & Dopita ( 2002), based on grids 
of pho toionization models, and the empirical one bv lBresohnl 
(2007), based on a sample of Hll regions with available 
r^-based abundances. 

(e) Ar3Q3 = log(f Ar ml A7 1 3 5/[Q hi] A5007), empirically 
calibrated by IStasins ka (2006). Although we measured the 
[Arm] A7135 line only in a minority of cases, we considered 
this indicator, monotonic as a function of O/H, because 
it does not involve the use of the [N II] A6583 line. G07 
proposed that the N/O ratio of H II regions in the extended 
disk of M83 could be higher than normal for their oxygen 
abundance, which would affect the oxygen abundances 
obtained from N2, 03N2 and [Nll]/[Oll]. While the results 
obtained from the auroral lines (Section [3. II ) do not support 
this idea, the Ar303 indicator could reveal obvious problems 
in abundances obtained from line ratios involving [Nil], and 
in the choice of upper/lower branch of R23. 

The results of the application of these strong-line indicators 
to our Hll region sample are shown in Fig.Hl where we plot 
O/H abundances as a function of galactocentric distance. In 
this figure w e also include comparison d ata for the inner disk 
of M83 f romlBresolin & Kennicuttl ( 120021 open light blue cir- 
cles) and lBresolin et al. (2005, light blue squares). The verti- 
cal dotted line in the diagrams represents the threshold radius 



Run (5.'1 = 6.7 kpc) defined bv lMartin & Kennicuttl d2001l) as 
the position where a strong break in the azimuthally averaged 
Ha surface brightness profile is encountered. Our first remark 
looking at the abundance trends in Fig. |4]is that the new Hll 
region data merge nicely with the comparison sample data. 
This is true also for the diagnostic diagrams in Fig. [5J thus 
excluding significant systematic offsets in our line fluxes. 

Qualitatively, the plots in Fig.|4]look similar to each other, 
although the scatter of the points varies for different indi- 
cators, being smaller in the case of [Nll]/[Oll]. The few 
datapoints available in the case of Ar303 (panel f) are in 
rough agreement with other indicators, suggesting that the 
choice of the upper branch of R23 (panel c), made on the 
basis of the strength of the [N II] A6583 line relative to Ha 
and [O II] A3727, is correct throughout the outer disk sample, 
and also confirms what is apparent from all the different in- 
dicators, i.e. that the oxygen abundances in the extended disk 
are not as low as 10%- 15% of the solar value, even though 
the absolute values really depend on the particular choice of 
abundance diagnostic. Concentrating now on the remaining 
diagrams, the scatter of the points is largest in the case of 
the N2 (panel a) and 03N2 (panel b) indicators. Looking at 
the inner disk data, both diagnostics seem to be affected by 
a satu ration effect at high abundance. This is a k nown ef- 
fect (IPettini & Pagell 12004 IKewlev & Ellison! 120081) . A log- 
linear upward trend in O/H towards the center of the galaxy 
is instead very clear in the case of the [N ll]/[0 II] plot. The 
familiar systematic offset between calibrations of abundance 
diagnostics based on photoionization model grids and em- 
pirical measurements is obvious in the comparison between 
panel d ([Nh]/[Q II] ind icator calibrated from the models by 
IKewlev & Dopital |2002[) and panel e (same indicator, cali- 
brated empirically by lBresohnll2007l) . In this case, the em- 
pirical method provides oxygen abundances that are nearly 
0.5 dex below the theoretical method. Higher abundances are 
also provided by the theoretical calibration of R23 (panel c) 
compared to the empirical calibrations. The radial trends are, 
however, virtually identical. 

The abundance offset between the N2 (panel a) and the em- 
pirical [N ll]/[0 II] (panel <?) diagnostics is puzzling at first, 
since both methods have been calibrated using samples of 
ionized nebulae for which the abundances were derived from 
the knowledge of the electron temperature (apart from a few 
objects at high metal licity used in the calibration of N2 by 
IPettini & Pagel 2004). To investigate the origin of this sys- 
tematic difference we ha ve considered the H II region sample 
used by Bresolin] (120071) to obtain his simple empirical rela- 
tion between O/H and th e r N Il"l/|"Q III line ratio, and verified 
that also in this case the IPettini & Pagell d2004l) N2 calibra- 
tion would provide abundances that are systematically ^0.2 
dex higher. T his might have to do w ith the inclusion of H II 
galaxies in the IPettini & P agel (2004) calibrating sample, and 
with the fact that a non-negligible diffuse [N II] -emitting com- 
ponent would affect the integrated spectra of galaxies in such 
a way as to lead to a slight o verestimate of the metallicity 
dMoustakas & Kennicutdl2006i) . On the other hand, the sim- 
ple analytical f orm of the empirical calibration provided by 
iBresolinl ((2007) could be slightly underestimating the O/H ra- 
tio in the range of interest here. Once again, however, we 
note that the general qualitative trends provided by different 
indicators are very similar. 

In panel e we have included the abundanc es d erived from 
the [Offl] A4363 auroral line, described in § [XT] for the four 
objects in Table [3] The good agreement between the direct 
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FIG. 4. — Radial abundance gradient in M83, where the oxygen abundance h as been determined from different strong-line me thods: (a) N2; ( b) 03N2; (c) 
R23; (d) [NlI]/[OlI] and the theoretical calibration by Kewley & Dopita I2002); (e) [NlI]/[OlI] and the empirical calibration by[Sresolin (200 7]); (f) Ar3Q3. 
The ne w M83 sample is represen ted by the red dots. The open circles and the full light blue squares represent the inner disk regions from Bresolin & K ennicutH 
1 2002) and Bresolin et al. (2005), respectively. The blue dotted line at R/R25 =0.79 shows the distance at which the azimuthally averaged Ha surface brightness 
has a sudden drop (Martin & Kennicutt 2001). In panel (e) we show with star symbols the four H II regions for which we have detected the [O III] A4363 auroral 
line (a line connects the 7" e -based abundances to the empirical ones). The linear regressions to the inner and the outer abundances (Eq. 1 and 2) are also shown 
with the dot-dashed lines in (e). 



abundances (star symbols) and the strong-line method abun- 
da nces indicates th at the empirical calibration of [N ll]/[0 II] 
by iBresoirnl (120071) is more tightly matched to the direct, T e - 
based method than the remaining diagnostics, at least for the 
objects considered here. 

In all the diagrams of Fig. [4] a break in the abundance gra- 
dient is clearly seen (except for Ar303, for which we do 
not have a sufficient number of points). Despite differences 
in the scatter in O/H and in the absolute value of the oxy- 
gen abundance, the gradient becomes virtually flat outside 
the isophotal radius, or slightly outside the threshold radius 
^Hn (blue dotted line). The effect is more clearly visible in 
the case of the abundances determined with the [N ll]/[0 II] 



line ratio, because of the smaller scatter. In fact, the dia- 
grams in panels d and e suggest the presence of two distinct 
groups of Hll regions, one located inwards of RJR25 = 1.2, 
for which there is a relatively steep decrease of O/H with ra- 
dius, and a second one from the isophotal radius outwards, 
for which the abundance is almost constant or very slowly 
declining with radius. In addition, a significant drop in metal- 
licity occurs around R/R25 = 1.2. Depending on the method, 
the oxygen abundance in the extended disk of M83 can be 
as high as 12 + log(0/H) = 8.6 (approximately solar, from R23 
and the theoretical calibration of [Nll]/[Oll]) or as low as 
12 + log(0/H) = 8.2 (1/3 solar, from the empirical [Nll]/[Oll] 
calibration). The auroral line-based abundances are in agree- 
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ment with the latter value. On the other hand, a solar-like 
metallicity at large galactocentric distances appears difficult 
to reconcile with the idea that the outer fringes of spiral disks 
are chemically unevolved. Finally, we note that a step-like 
gradient in the outer disk of M83 is suggested by the abun- 
dances published by G07 (their Fig. 7a), although the abso- 
lute abundances are smaller than ours, and virtually only one 
out of a total of 12 Hll regions with a published metallicity 
(XUV 01) is located well beyond the isophotal radius, with an 
oxygen abundance 12 + log(0/H) = 7.86. 

In order to quantify the abundance gradient, we have carried 
out a linear reg ression to the data points shown in Fig. |4je) 
dBresolinll2007l calibration of [N ll]/[0 II]), independently for 
the inner and the outer disk. We have arbitrarily included in 
the inner disk fit the five H II regions that lie at R > R25 and 
12 + log(0/H) > 8.3 (all relatively close to each other south 
and south-west of the center), since also visually they fit very 
well an extension of the inner gradient. We obtained: 

• inner disk (54 data points, rms = 0.05): 

12 + log(0/H) = 8.77(±0.01)-0.25(±0.02)/?//? 25 (1) 

• outer disk (19 data points, rms = 0.05): 

12 + log(0/H) = 8.25(±0.04)-0.04(±0.02)/?//? 25 (2) 

The same slopes, but with a systematic offset of +0.47 
dex, are obtained using t he [N ll]/[0 II] abundances from the 
iKewle v & Dopita (2002) calibration. The two regressions in 
(1) and (2) cross at R/R25 = 2.26. The slope of the gradients 
correspond to -0.030 dexkpc" 1 (inner) and -0.005 dexkpc" 1 
(outer). The inner abundance gradient derived here is slightly 
steeper than the value reported earlier (-0.02 2 dexkpc" 1 : 
iBresolin & KennicuttH2002l: iPilvugin et al.ll2006l) . The slopes 
obtained by G07 are very different, due to the discrepancies 
in the comparison with our abundance analysis (§ I3.41 >. and 
to the assumed shape of the gradient. Expressing the inner 
gradient slope in terms of the isophotal radius, our value of 
-0.27 dex//?25 is sim ilar to the mean of - 0.23 dex//?25 for 
five barred galaxies in Zaritsk v~et al.l (119941) . compared to the 
mean for their whole sample (barred + unbarred spirals) of 
-0.59 dexkpc" 1 . We stress, however, the difficulty of mea- 
suring chemical abundance gradients in galaxies, and espe- 
cially the fact that the use in the literature of different metal- 
licity diagnostics leads to very different results. An example 
is provided by the nearby galaxy M33, for which the estimates 
of the slope of the oxygen abundance g radient vary between 
-0.012 dexkpc" 1 (ICrockett et al.ll2006l) and -0.10 d ex kpc" 1 
dVflchez et aklll988l see lRosolowskv & Simonll2008l) . 

In order to better visualize the spatial distribution of the 
oxygen abundance, and in particular in relation with the lo- 
cations of the UV-emitting knots, we show in Fig. [5] and [6] 
the positions of the H II regions in the northern and southern 
fields observed by us on top of the GALEX image of M83 4 . 
We labeled each H II region (represented by squares) with the 
12 + log(Q/H) valu e obtained from the [Nll]/[Oll] ratio and 
the Bresolin (20 071) calibration , as well as from the N2 index 
and the Pettini & Pagel (2004) calibration (values in brack- 
ets). In Fig. |5]the sudden drop in oxygen abundance that oc- 
curs beyond the optical isophotal radius (the dashed curve) is 
impressively displayed. The six H II regions at the inner edge 
of the optical disk have an average 12 + log(0/H) = 8.59 (8.75 
if N2 is used), with a very small dispersion, while just beyond 

4 Available at http://photojournal.jpl.nasa.gov/catalog/PIA10374 




FIG. 5. — The positions of the northern field H II regions shown as squares 
on the GALEX image of M83 (courtesy NASA/JPL-Caltech), an d labeled 
with t heir 12 + log(0/H) values obtained from [NlI]/ [OlI] in the [Bresolin 
120071) scale, and from N2 in the Pettini & Pagel ( 2004) calibration (in brack- 
ets). The nucleus of the galaxy appears at the bottom. The dashed curve 
represents the R25 radius. 

this radius the value drops by ~0. 30-0. 35 dex, and remains 
virtually constant along the filament that reaches the bright 
cluster of UV knots at the top of the image. Qualitatively the 
same picture is obtained if an alternative metallicity diagnos- 
tic, such as R23, is used. It is also worth pointing out that the 
two groups of H II regions on opposite sides of the R25 bound- 
ary line have similar H/3 luminosities and [O lll]/[0 II] ratios, 
which implies that the drop in O/H across the optical edge 
cannot be attributed to systematic differences in ionizing flux 
output or nebular excitation. 

In the southern field (Fig. |6]) the situation could be more 
complex, depending on the choice of empirical abundances. 
While also in this case the inner H II regions near the edge of 
the bright optical disk have a mean abundance that matches 
very well that obtained at the opposite side of the disk 
[12 + log(0/H) = 8.62 using [N ll]/[0 II], notice the small vari- 
ation from this mean for the 10 objects in Fig. [6), outside R25 
there appears to be a number of nebulae at an intermediate 
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FIG. 6. — The positions of the southern field H II regions shown as squares 
on the GALEX image of M83 (courtesy NASA/JPL-Caltech), an d labeled 
with t heir 12 + log(0/H) values obtained from [N II]/ [Q II] in the [Bresolin 
120071) scale, and from N2 in the Pettini & Pagel ( 2004) calibration (in brack- 
ets). The dashed curve represents the R25 radius. 

oxygen abundance between that of the inner disk and the outer 
disk, with 12 + log(0/H) ~ 8.4. These belong to the group of 
points in Fig. |4] that lie along the exponential fit to the inner 
gradient. Within the same large clustering of UV knots (left 
of the bright foreground star) some H II regions have oxygen 
abundances that are 0.2 dex smaller, matching those found at 
the end of the filament near the bottom of the image. The N2 
abundances (in brackets in Fig. [6]l provide, however, a more 
uniform set of abundances beyond /?25- A possibile explana- 
tion for the observed variations can be the presence of aper- 
ture effects affecting some of the objects, since for this field 
the red and blue spectra were obtained on separate nights. 

3.3. Nitrogen 

In Fig. [7] we show the radial trend of the N/O ratio (top) 
and the N/O ratio as a function of O/H (bottom). N/O 
has been measured following the algorithm proposed by 
iThurston et al.l d 19961) . which provides the electron tempera- 



fa = 0.6065 + 0. 1600x + 0. 1878x 2 + 0.2803x 3 (3) 



(x = logR23, t2 in units of 10 4 K), and equation (9) by 
iPagel et ail (1 19921) used to derive N + /0 + from the observed 
[Nll]/[Oll] line ratio: 



, N + , [Nil] AA6548. 6583 

log — = log 

5 0+ 5 [OII]A3727 



-0.307-0.021og? 2 - 



0.726 



(4) 



We then assume N + /0 + = N/0 As the top panel of Fig. [7] 
shows, the N/O ratio is virtually flat around N/O ~ -1.2 
(~0.5x the solar ratio) outwards of the isophotal radius. 
From the bottom panel of Fig. [7] where the O/H abun- 
dance is obtained from the N2 diagnostic (to avoid using 
the [N ll]/[0 II] ratio for both the O/H and N/O abundance 
ratios), we infer that the N/O ratio as a function of O/H 
throughout the disk of M83 is consistent with the obser- 
vations in other nearby galaxies, such as the well-studied 
M 101 (squares, corresp onding to the r ^-based abundances 
of iKennicutt eTaIll2003l) and NGC 2403 dGarnettet al.ll 19971 
green triangles). The N/O abundance ratios based on the 
[O III] A4363 auroral-line (open star symbols) are in good 
agreement, within the scatter, with those determined from 
R23 (also in the case of the M101 Hll regions the agree- 
ment is on the order of 0.1 dex). The N/O ratio is approx- 
imately constant below 12 + log(0/H) = 8.4, then it rises lin- 
early with O/H. Qualitatively, this is explained in terms of a 
primary production of nitrogen, which gives rise to the con- 
stant pedestal in N/O at low oxygen abundance, and of a sec- 
ondary co mponent, which is proportional t o the oxygen abun- 
dance (e.g. IVila Costas & Edmunds 1993). We note that the 
oxygen abundance scale obviously depends on the choice of 
metallicity diagnostic. 

3.4. Comparison with previous work 

The work by G07 was the first to characterize the chemical 
abundances in the XUV disk of M83 (as well as NGC 4625). 
These authors reported line emission fluxes for 22 H II regions 
in this galaxy, 19 of which lie in the outer disk. The chemical 
analysis was based on a subset of 12 objects (9 in the outer 
disk), for which reliable fluxes in both the [O II] and [O III] 
lines could be obtained, therefore allowing the determination 
of the R23 parameter. Lacking [O III] A4363 detections, G07 
estimated the oxygen abu ndan c es fro m the R23 method, using 
the ca librations of iMcGaughl d!991l) and iPilvugin & Thuaij 
(2005). Apart form the three innermost H II regions, G07 pro- 
vide both upper and lower branch abundances for all objects, 
since, according to these authors, the R23 degeneracy was dif- 
ficult to remove. While the [Nll]/[Oll] and the [Nll]/[Sll] 
line ratios would suggest abundances between Z Q /2 and Z Q 
(consistent with upper branch abundances), their preferred 
abundances are between Z Q /10 and Z Q /4 (lower branch), 
based on photoionization modeling in which the N/O ratio 
is constant as a function of O/H and equal to the solar value 
(models with an O/H-scaled value for the N/O ratio yielded 
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FIG. 7. — (Top) Galactocentric gradient of the N/O ratio, expressed 
in terms of the isophotal radius Rjs- Symbols as in Fig. [4] Lines 
connect empirical abundances with T e -based abundances. (Bottom) 
Relationship between the N/O ratio and the O/H ratio. The latter has 
been determi ned via the N2 diagnostic. A comparison sa mple has been 
drawn from Kennicutt et al. 12003, M101, squares) and Garnett et al. 
119971 triangles). The solar symbol i s draw n at 12 + log(O/H)0 = 8.66 
and log(N/O) =-0.88 lAsplu nd et al. 20 0$). Fo r M 10 1 we plot the 
original T e -based abundances from Kennicutt et al. 12003, full squares) 
and those derived from N2 (for O/H) and Eq. 4 (for N/O, small open squares). 



less satisfactory fits to the emission line ratios). This produces 
a discrepancy in relation to our work: instead of measuring 
a mean abundance in the outer disk of 12 + log(Q /H) rj 8.0 
(Z©/5; or lower if using the Pilyugin & Thuan 2005 method), 
we measure abundances that are between 0.2 dex and 0.6 
dex higher (Zq/3-Zq), depending on the metallicity indica- 
tor used. As a cautionary note, we point out that the O/H 
abundance ratios derived from N2, 03N2 and [N ll]/[0 II] are 
calibrated using H II region da ta and models with a 'standard ' 
dependence of N/O on O/H (IVila Costas & Edmundslll993l) . 
A higher N/O ratio in the outer disk, as proposed by G07, 
would result in smaller O/H values. Our preferred abundances 
in the outer disk are those that are tied to the [O III] A4363 re- 
sults, i.e. 12 + log(0/H) «8.2. Our results on the N/O ratio 
(also those based on the [O III] A4363 detections) do not sub- 
stantiate the choice made by G07 of a solar N/O ratio in the 
outer disk of M83: the derived N/O ratios appear normal for 
the given O/H ratio, given the scatter observed in extragalac- 
tic H II regions (the 7^-based N/O ratios appear slightly higher 
than the average). Our results are largely independent of the 
R23 degeneracy, since we have used metallicity indicators that 
are mono tonic with oxygen abundance. 

4. DISCUSSION 

The main findings of our observations of the extended disk 
of M83 can be summarized as follows: (a) a flat oxygen 
abundance gradient beyond R25, with a relatively high O/H 
ratio; and (b) a discontinuity in the oxygen abundance ra- 
dial distribution at the edge of the optical disk. The latter 
is most obvious when the oxygen abundance is derived from 



the [Nll]/[Oll] ratio. The discontinuity takes place where 
the radial abundance gradient turns from exponential to fiat. 
This radical change in the abundance distribution just beyond 
the optical radius of M83 implies different star formation his- 
tories between the inner and outer parts of the disk of this 
galaxy. 

These findings are supported by new emission line data for 
nearly 50 H II regions, distributed in galactocentric distance 
between 0.6 and 2.6 isophotal radii, and by the application of 
different metallicity diagnostics, including the [O III] A4363 
auroral line and several strong-line methods. In this Section 
we first look at previous detections of changes in the slope 
of the abundance gradients (what we otherwise call breaks) 
in other galaxies. We then critically analyze our abundance 
data, in order to test whether our results depend on the choice 
of metallicity diagnostics. Finally, we attempt to provide an 
interpretation for the observed features in the abundance dis- 
tribution. 

4.1. Bimodal abundance gradients: observations and theory 

Changes in galactic abundance gradients, with shallower 
(or flat) slopes at large galactocentric distances, have already 
been determined from studies of the H II regions located in 
the main optical disk of a few spiral galaxies. For example, 
iMartin & Rovl £l995) found a break in the slope of the abun- 
dance gradient of the barred galaxy NGC 3359, in correspon- 
de nce of the corotation radius. A similar break was observed 
bv lRov & Walshl d 1 9971) in NGC 1365, also a barred galaxy. In 
both these cases, the break in the gradient occurs well within 
the bright optical disk of the galaxy. In M83, instead, we find 
that the flattening occurs just outside of the main star-forming 
disk, approximately at the isophotal radius. 

A bimodal chemical abundance distribution is observed 
also in the Milky Way. The Cepheid abundances of vari- 
ous elements, and iron in particular, show a discontinuity (on 
the order of 0.2 dex) in their radial distribution, with a shal- 
lower gradient at galactocentric distances larger than 10 kp c 
dLuck et alj|2003t lAndrievskv et al.1 12001: lYong et al.|[2 006). 
A sim ilar discontinuity has been observed by iTwarog et alj 
d 1997b in the metallicity distribution of Galactic open clus- 
ters. The lack of efficient radial mixing beyond the coro- 
tation resonance (where the spiral pattern speed equals the 
galactic rotation velocity) and the gap in gas density at 
corotation could be responsible for the break in the abun- 
dance gradient and for t he discontinuity determined from 
the C e pheid observations (|Mishurov et a l. 2002; Lepi ne et alj 
120031: lAcharova et alj l2005h . Th e behavior of the ionized 
gas abundances is less clear, since iVflchez & Estebanl (fl996l) 
obtained an essentially flat gradient in the nitrogen and 
oxygen abundanc e s bey ond 12 kpc from the center, and 
Maciel & Ouireza (1999) measured a similar flattening for 
oxygen and other elements from plane tary nebulae located be- 
yond 10 kpc, while o ther authors (e.g. Deharveng et alj[2000b 
iRudolph et alj |2006) found no evidence for breaks in the H II 
region abundance gradient. 

Turning to the theoretical side, non-linear gradients (in log- 
arithmic abundances), flatter at large galactocentric distances, 
can be reproduced by 'inside-out' scenarios of galaxy for- 
mation in which the Milky Way disk is built up via gas in- 
fall, and in which the timescale for the formation of the disk 
increases with galactocentric dist ance (Matteucci & F rancois! 
[19891: iBoissier & Prantzosl fl999l; IChiappini et alj|2003l) . In 
these models the radial variations of the gas infall rate and 
of the star formation rate in the galactic disk control the for- 
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m ation and the slope of t he abundance gradients. As shown 
by lChiappini et alJ ([2001), the shape of the abundance gradi- 
ents at large galactocentric distances is regulated by the evolu- 
tion of the halo. Despite the fact that different authors predict 
that the overall disk abundance gra die nts become eithe r flat- 
ter (e.g. [Prantzos & Boissier 2000 and lHou et al.ll206ol who 
also predict a flatter gradi ent in the inner disk) or steeper 
(e.g. IChiappini et al.l 1200 lb with time, recent galactic evo- 
lution models are able to provide a good match to the ra- 
dial trends of many different chemical species measured fo r 
Cepheids, young stars and open clusters (ICescutti et al.l20O7b . 
The development of a plateau in the abundance gradient at 
large galactocentric distance (as well as in the central re- 
gion) is also predict ed by the chemodynamical models of 
Samla ncTet al.l (119971) . For galaxies beyond the Milky Way, 
however, it is hard to generalize these model results, since flat 
abundance gradients and/or discontinuities are rarely, if ever, 
observed. In fact, almost no information on nebular abun- 
dances beyond the optical edges of spiral disks is available. 
In summary, changes in the abundance gradients, with flat- 
ter slopes in the outer regions, as well as abundance discon- 
tinuities, have been observed in other galaxies, although very 
rarely (a discontinuity only in the Milky Way). Despite some 
qualitative similarities with the examples mentioned above, 
we cannot conclude that the physical causes of the abundance 
distribution in M83 are necessarily the same. 

4.2. Is the abundance break in M83 real? 

iPilvuginl (120031) has questioned the reality of slope changes 
in the radial abundance gradients of spiral galaxies, showing 
that an artificial break can occur when abundances are deter- 
mined via the R23 method, either as a result of adopting the 
wrong branch of the indicator outside of a certain galactocen- 
tric distance, or as due to a break in the radial trend of the 
excitation parameter P = [O m]/([0 ll]+[0 III]). We think that 
our result in the case of M83 is quite robust, having inferred 
the abundance break from a variety of metallicity indicators, 
most of which are monotonic with O/H, thus being unaffected 
by the degeneracy of the R23 method. To strengthen our con- 
cl usion, we procee ded with a test similar to that carried out 
by Pilvugin] (120031) . i.e. we applied our abundance i ndicators 
to the Hll regions of the galaxy M101, for which Pilvugin 
d2003l) showed that the use of the R23 indicator can lead to 
a fictitious break in the gradient. The results of our test are 
shown in Fig. [8] In panel (a) the oxygen abundance has been 
obtained from the R23 diagnostic a nd the P-method, as cali- 
brated bv lPilvugin & Thuanl (12005). The green symbols refer 
to the upper branch calibration, while the light blue symbols 
refer to the lower branch calibration. Full dots are used for 
N2 >— 1, open circles for N2 < — 1, This plot shows that, in- 
deed, the adoption of a single, upper branch calibration of R23 
throughout the whole disk produces a radial trend that devi- 
ates in the outer galactic region from the linear gradient mea- 
sured from the [O Hll A4363 aur oral line (triangles and dashed 
line, from lKennicutt et aill2003l) . here chosen to represent the 
real radial trend of the oxygen abundance. The break occurs 
approximately at the radius in which a break in the P parame- 
ter occurs [panel (d)]. The N2-based criterion to discriminate 
between the upper and lower branches appears effective at the 
extremes of the abundance range, while in the intermediate 
range 8.0 < 12 + log(0/H) < 8.3 the branch selection from N2 
is less in agreement with the A4363-based results. This is not 
unexp ected, since the P-metho d is not calibrated within this 
range (Pilvu gin & Thuanll2005l) . 



The abundances in panel (b) of Fig. [8] were obtained 
from t he [N n]/[0 11] ratio, a s calibrated bvlKewlev & Dopital 
(120021 light blue circles) and lBresolir] (120071 red dots). In the 
latter case, we used the larger symbo ls for H II regions wit h 
available [Olll]A4363 abundances dKennicutt et all 120031) . 
Apart from the systematic offset noted previously, the two cal- 
ibrations provide virtually the same abundance gradient, with 
a slope that is in good agreement with the A4363-based one, 
and without any obvious deviation from a linear trend. A pos- 
sible decrease in the gradient slope of the strong-line abun- 
dances, suggested by the single H II region at R/R25 = 1 -25 , 
might occur at 12 + log(0/H) < 8.0 in the iBresolirj (120071) 
scale, i.e. well below the abundance at which the gradient in 
M83 levels off to a constant O/H. This suggests that the rela- 
tive abundances determined in M83 from the [N ll]/[0 II] ratio 
are reliable throughout the inner and the outer disk. The same 
conclusion can be reached from panel (c) of Fig. [8] where 
the N2 indicator wa s used. We sh o w the abundances from 
the calibration s of [P ettini & Pagej (120041 green dots) and 
iDenicolo et all d2002l open light blue circles). Once again, 
the agreement with the abundances from the auroral-lines is 
good down to 12 + log(0/H) = 8.2, while the N2 abundances 
in the outer disk of M83 flatten around 12 + log(0/H) = 8.4 
(see Fig. |4). In conclusion, this test supports our results con- 
cerning the abundance gradient in M83, and in particular the 
fact that the slope change in the abundance gradient seen at 
the isophotal radius is real. 

4.3. Explaining the abundance distribution in M83 

If the abundance gradient in M83 has the shape shown 
in Fig. |4] with both a flat trend and a discontinuity occur- 
ring around the isophotal radius, it is important to under- 
stand how these features can be explained within a model 
of galactic evolution. The fact that nowhere in the extended 
disk of the galaxy the oxygen abundance is below 1/3 so- 
lar (or possibly even higher, depending on the abundance 
diagnostic that is used) is also remarkable. Because of the 
low levels of star formation activity at large galactocentric 
distances, as measured from the averaged Ha profile, we 
would expect low rates of chemical enrichment and there- 
fore rather low oxygen abundances in the outskirts of the 
galaxy. Our findings go against this simple expectation, in- 
dicating that some processes are at work that both homog- 
enize and increase the chemical abundances even far from 
the region of main galactic star-forming activity. However, 
a non-negligible level of chemical enrichment in extended 
spiral disks is expected, based on the observed UV emis- 
sion (which measures star formation on timescales on the 
order of a few hundred Myr), as well as on the presence 
of older generations of stars that are resolved beyond the 
optical edges o f nearby late-type spirals (e.g. NGC 2403 , 
NGC 300, M33: lDavidgdl2003l iBland-Hawthorn et aT]|2005l; 
iBarker et al.l 120071). Moreover, the discovery of inter galac- 
tic H 11 regions dGerhard et al.ll2002t iRvan- Weber et aU2004l) 
has led to the determination of near-solar gas-phase oxy- 
gen abundances at pr ojected distances up to 25-30 kp c from 
the nearest galaxies (M endes de Oliveira et al.ll2004l) . In at 
least some of these cases the star forming activity appears 
to have been triggered by tidal interactions between galaxies, 
which could also help explain the high metallicity as a result 
of chemicall y enriched gas be ing stripped out of metal-rich 
galaxy disks (ISakai et al.ll2002l) . A crude estimate of the ex- 
pected present-day gas abundances in the outer disk of M83 
can be obtained by assuming that star formation has been 
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FIG. 8. — Abundance gradient in M101 from different strong-line indi- 
cators: (a) R23, showing the res ults of the upper branch (green) and lower 
branch (light blue) calibrations of Pilyugin & Thuan (2005). Full dots corre- 
spond to values of t he N2 indicator larger than —1 .0. (b) The [N I1]/[ Q II] indi- 
cator, calibrated by Kewlev & Dopita 12002, light blue circles) and Bresolin 
(2007, red dots). Larger symbols represent Hll regions with available 
[O III] A4363 auroral line measu rements, (c) N2, cal ibrated bv lDenicolo et al.l 
12002, light blue circles) and Pettini & Pagel 12004, green dots). Larger sym- 
bols represent H II regions with available [O III] A4363 auroral line measure- 
ments. The triangles represent H II regions with oxygen abundances obtained 
from the A4363 line. The dashed line is the abundance gradient obtained 
from a linear fit to these abundances, (d) T he radial trend o f the P parameter 
in M1 1. The observationa l dat a are from Kennicutt et al. (2003), Bresolin 
<2007T) ,rvan Zee et all (19981) and lKennicutt & Garnett 1 1 99(3) 1 



ongoing for at least the last Gyr, as deduce d by iDong et al.l 
(2008). However, numerical simulations bv lBush et al.l (2008) 
suggest that star formation is a long-lasting process in XUV 
disks, extending over timescales that can be as long as the life- 
time of the galaxy. We have estimated the current SFR density 
(Esf r) of the outer disk of M83 by measuring the far-UV flux 
detected in the form of UV clusters and complexes beyond 
R25 and out to 4/?25- The value we found, after correcting 
for the presence of older clusters and background sources, is 
^sfr = 1-0 x 1O~ 5 M yr" 1 kpc -2 , much l ower than the esti- 
mate of 8 x 1 0" 4 M Q yr" 1 kpc" 2 obtained bv lDong etail (120081) 
from the 8 /im emission in two Spitzer IRAC fields at approx- 
imately 2.3 R25 from the center. We attribute much of the dis- 



crepancy to the fact that we consider the SFR integrated over 
the whole outer disk (out to 4 ^25), therefore averaging the 
far-UV emission from star-forming complexes with the null 
contribution from large zones where no clusters are detected. 
The two fields in Dong et al. (2008) were selected based on 
the presence of active star forming regions, and have there- 
fore an enhanced Y,sfr relative to the mean of the outer disk. 
Our lower value is thus more representative of the SFR of the 
outer disk as a whole, and would be more appropriate to adopt 
when considering an evolutionary scenario in which sites of 
star formation appear stochastically throughout the extended 
disk, perhaps a s a consequence of triggering by outer spiral 
perturbations (Bus h et 



sequence 
SE008). 



With the current rate of star formation, as determine d 
above, and a net oxygen yield of 0.01 ( Ma ederl |1992|) . 
~10 2 M Q of oxygen would have been released per kpc 2 
into the interstellar medium of the outer disk over the past 
Gyr. With a gas surface density E(H I) ~ 1.6M w pc~ 2 at 3 #25 
(derived combining informat i on f rom Rogstad et al. 119741 
i Huchtmeier & Bohnenstengefl 119811 and ICrosthwai te et al.l 
|2002|) . we then estimate an O/H ratio roughly equal to 1/7 of 
the solar value [12 + log(0/H) = 7.8] as a result of the chem- 
ical enrichment over the past Gyr. Considering the approx- 
imations and uncertainties involved in this estimate, we can 
therefore reasonably recover the oxygen abundance measured 
in the extended disk of M83 by assuming a constant SFR over 
the past 2-3 Gyr. Longer timescales, up to several Gyr, can be 
easily accomodated, for example by assuming lower SFR in 
the past than the cur rent value, or smaller oxygen yields (as in 
Pilyugin et al. 2007). This demonstrates that a low, but persis- 
tent, level of star formation can lead to a substantial chemical 
enrichment, compatible with the observed abundances. We 
stress that this association between gas density, star formation 
and metallicity can only be made in regions where star for- 
mation has been taking place. There may well be areas of the 
outer disk that have sub-critical gas densities and are there- 
fore quiescent and metal poor if they have remained below 
the critical density during most of their history. 

We also point out that, considering the abundances of 
the outer disk of M83 within the context of a closed box 
chemical evoluti on model, the effective yield y e ff = Z/(ln /i -1 ) 
(Edmunds 1990) would be quite large in the extended disk of 
M83, considering that the gas mass fraction [i must be near 
unity (although no estimate of the stellar mass is available, 
its fractional component is presumably very small, given 
the absence of an easily detectable stellar disk). In other 
words, the outskirts of M83 are overabundant for their gas 
fracti on, the opposite of what is observed in dwarf gal axies 
(e.g. lMatteucci & Chiosill 19831: Ivan Zee & Havnesll2006l) . 

As noted above, gas-phase abundance plateaus have been 
inferred for a few barr ed galaxies from the measurement 
of stro ng-line indexes dMartin & Rovl 119951; iRov & Walshl 
119971) . In these cases the explanation of the observed 
breaks in abundance gradients is given in terms of the non- 
axisymmetric potential of a young and strong bar, in which 
radial gas flows produce, with time, increasingly shallower 
slopes across the galactic disk, as a consequence of the re- 
distribution of mass and angular momentum, and a 'steep- 
shallow' transition at the galactocentric distance that corre- 
sponds appro ximately to the corotation radius (Fried iiet al.1 
ll994t [Friedh' 1998). Gas flows that are preferentially directed 
outwards beyond corotation have a diluting effect, producing 
a flattened abundance distribution in the outer part of the disk. 
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The bar in M8 3, with a deprojected axis ratio (b/a) = 0.38 
(Martin 1995), is among the strong est (lowest b Id) in the sam - 
ple of barred galaxies studied by Mar tinet & Friedlil d 1997b . 
who show that stronger bars corre late with shallower r adial 
abundance gradients. According to Mar tin & Royl (1 19951) . the 
weak star formation in the bar, the strong circumnuclear ac- 
tivity and the shallow abundance gradient, as observed in the 
disk of M83, support the idea that this is a rather evolved 
barred system. Within a picture of secular evolution of barred 
galaxies, an initial phase of intense star formation, that pro- 
duces a discontinuity in the abundance gradient at the edge 
of the bar, is followed by an almost quiescent phase in the 
bar, during which the funneling of gas toward s the central re - 
gions trigge rs star formation near t he nucleus (Combes 2008). 
Recent CO dLundgren et alj|2004l) and Ha dFathi et alj|2008l) 
kinematical data have confirmed the presence of a gas inflow 
along spiral trajectories, from kpc-size scales down to a few 
tens of pc from the nucleus of M83, and which is likely driven 
by the potential of the bar. 

While the gas flows induced by the presence of the bar 
are known to generate overall flatter abundance gradients 
in the in ner disk of spiral galaxies compared to non-barred 
galaxies dVila-Costas & Edmunds|[l992l: IZaritskv et alj|1994t 
lDutil&Rovl fl999). this mechanism seems unlikely, or at 
least insufficient, to explain the abundance trend in the outer 
disk of M83. As mentioned above, the abundance break 
caused by a bar is predicted to occur at the galactocentric 
distance that corresponds to corotation. In the case of 
M83, corotation is at 2.8 3 arcmin from the center (3.7 
kpc, iLundgren et alJ 120041) . while the break in abundance 
gradient occurs much further out, at the edge of the optical 
disk (7?25=6.44 arcmin, 8.4 kpc). On the other hand, 
angular momentum redistribution and gas fl ows predicted 
for v i scous evolution of star-forming disks dLin & P ringle 
fl987t lOarkel [19891 iFerguson & Clarkd l200lh can give rise 
to fla t gradients in the ISM of the outer parts of spiral 
disks dSommer-Larsen & Yoshiilll990tlTsujimoto et alJl!995b 
iThon & Meusingerlll998lT ~ 

An alternative explanation can be explored, by noting that 
in M83 the drop in the azimuthally averaged Ha surface 
brightness at Z?Hn occurs in proximity of t he outer Lind- 
blad resonance, at 4.83 arcmin (6.3 kpc, Lundgren et al. 
2004), which could act as a possible potential barrier for 
gas flows, and wh ich might cause discon tinuities in the 
star formation rate dAndrievskv et alJ 120041) . We also note 
that the abundance trends shown in Fig. [4] mimic the to- 
tal gas (neutral + molecu lar) surface density profile shown by 
Crosthwait e et al.l ([2002, their Fig. 8), which flattens at ap- 
proximately 6 arcmin from the center, beyond the sharp de- 
cline in both CO and H I emission at the edge of the optical 
disk (but we must point out that molecular data have not yet 
been p ublished for the disk outside /?2s)- ICrosthwaite et alJ 
(2002) pointed out that the abrupt discontinuity between the 
molecular-dominated inner disk (r < 5') and the atomic outer 
envelope is rather unique among late-type spiral galaxies. In 
relation to this finding, a modest drop in the mean hydrogen 
volume densities of giant molecular clouds a t R > R25 has re- 
cently been measured by iHeiner et al.l d2008l) . Therefore, the 
flat abundance gradient we observe in the outskirts of M83 
may simply be a reflection of the gas density profile beyond 
the optical radius. This is illustrated in Fig. [9] where we show 
the empirical abunda nces taken from Fig. CT e), and the gas 
density profiles from Rogsta det al.l d 1974 their model Hi 




R/R 25 

FIG. 9. — Comparison between the empirical O/H gradient (from Fig.|4j) 
and the logarithmic gas density profile (units for the latter are shown at the 
top right). The two to p curves represent th e surface density distri bution of 
H i (f rom the model of lRogstad etalJI 19741) and Ho (from Crosth waite et all 
12003) . 

distribution in Fig. 7) and Crosth waite et alJ d2002l their H2 
data). As the comparison shows, the relative slopes of the 
inner and outer O/H distributions closely match the relative 
slopes of the £(H2) and E(H I) profiles, respectively. The dis- 
continuity in the O/H distribution occurs where the H I surface 
density starts to dominate the radial gas profile. 

Beyond the optical edge, a large envel o pe of 
neutral hydrogen is pre s ent dRogstad et alJ Il974t 
iHuchtm eier & Bohn enstengell 119811) . extending 47 ar- 
cmin (62 kpc) from the c enter, and containing an Hi ring 
outside the optical radius (Tilanus & Allen 1993). The low 
gas surface density in the envelope agrees with the low 
level of star formation in the outer disk. Under such a 
situation, with low efficiencies of the processes that lead to 
star formation, flat abundance gradients are expected, due to 
the unevolved statu s of the system (a s shown, or example, 
by the models of iMolla & Dfazll2005l) . The situation might 
be similar to that found in low surface brightness galaxies 
(LSBGs). A possible analogy between XUV disks and the 
gas-ric h LSBGs has already been noted by Thilk eTet alJ 
(2007). The star forming activity in LS BGs takes place at 
low, p ossibly subcritical, gas densities dvan der Hulst et ail 
119931: Ide Blok et ail 1 1996b . Low rates of star formation, 
perhaps occuring sporadically or inter mittently, have been 
obtained for LSBGs by several authors dGerritsen & de Blokl 
119991: Ivan den Hoek et all 120001: iBoissier et al.l 120081) The 
low chemical abundances of LSBGs dMcGaugh 1 19941) 
denote fairly unevolved systems, even if there are exceptions 
(Bergmann et alJ 120031). The notion of LSB Gs as 'slowly 
evolving' svstems lMcGaugh & de Blokll 1997b might then fit 
the description of outer spiral disks, as we ll. The H II regio n 
oxygen abundances meas u red in LSBGs by McGaugh ( 1994), 
de Blo k & van der Hulstl d!998l) and iKuzio de Narav et alJ 
(120041) span a wide range. 12 + log(0/H) = 7.6-8.8, without 
a preferred value (these O/H abundances, obtained from R23, 
are in the lMcGaughl ll99ll theoretical scale). The oxygen 
abundance value obtained in the extended disk of M83 lies 
comfortably within this range and, as shown earlier, is con- 
sistent with the enrichment level one obtains with the present 
SFR in the outer disk over timescales of a few Gyr. Moreover, 
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the work by de Blo k & van der Hulstl 0998) showed that the 
radial distribution within LSBGs galaxies is approximately 
flat, i.e. there is no evidence for the presence of gradients in 
the oxygen abundance. As suggested by these authors, this 
finding would agree with the sporadic star formation rates 
measured in LSBGs. Star formation in low surface brightness 
disk galaxies, as well as in the outer disk of M83, seems to 
favor a chemical evolution in which large-scale gradients are 
not produced. 

The wide variety of metal enrichment seen in LSBGs, as 
well as dwarf galaxies, may reflect their range of morpholo- 
gies: some ha ve a quite obvious e xtended disk structure 
(e.g. NGC6822, Weldrake et ai1l200l . while others are much 
more irregular. Perhaps it is the LSBGs with rotationally 
supported disks that can maintain slow, but steady star 
forma tion. In the case of the outer disk of M83, iDong et all 
(2008) find that the Hi surface density distribution peaks at 
the critical density (but we note again that the total H2 + H I 
surface density of the outer disk is not known, due to the lack 
of molecular data). 

Finally, we should also entertain the idea that peculiari- 
ties in the chemical abundance distribution of the extended 
disk of M83 could arise from the gravitational interaction 
with one or more o f the several dwarf gala xies found in 
the vicinity of M83 ( Karache ntsev et al]|2007l) . An interac- 
tion with NGC 5253 about 1 Gyr ago has been p roposed by 
Ivan den Berghl (Il980i see alsolCalzetti et al.lll999l). T he faint 
dwarf KK 208 dKarachentseva & Karachentsev! 1998h . 25 kpc 
in projection to the north of the center of M83, is strongly 
elongated, likely as a resu lt of the tidal interaction with M83 
( Kara chentsev et al]|2002l) . Its location does not correspond 
to any of the UV-bright filaments discovered by GALEX, but 
UV knots are visible at even larger galactocentric distances. 

The warped structu re of the extended H I envelope of M83 
(Rogst adet al.lfT974l) is a possible signature of past galaxy 
encounters. Moreover, from the kinematic point of view 
the outer Hi envelope appears to be decoupled from the 
inner part, with disk rotation taking place on different planes 
(Huchtmeier & Bohnenstengel 1981). Gas stripping and/or 
redistribution taking place as an effect of the interactions 
could help explain the different chemical properties betweeen 
the inner and the outer part of the disk. For instance, we could 
explain the high metal abundance of the outer disk if we 
postulated that it is mainly composed of gas that was tidally 
transported from chemically evolved regions at smaller radii 
in the main disk of M83. 

5. CONCLUSIONS 

The outskirts of spiral galaxies hold a powerful diagnostic 
value in the study of their formation and evolution. Breaks 
in the exponential surfac e brightness profiles are detected up 
to six disk scale lengths (Pohlen & Trujillo 2006; Erwi n et all 
2008), reflecting at least in some cases breaks in the mass 



distribution dBakos et alj|2008l) . These data are spawning a 
number of models for their interpretation, and include the 
effect of radial star formation drop s resulting from drops in 
the gas density dRoskar et al.ll2008l) . star formation thresholds 
(lElmegreen & Hunte r 2006) a nd angular mome ntum redistri- 
bution by galactic bars dDebattista et al.ll2006l) . In the UV- 
detected extended spiral disks star formation takes place even 
beyond these breaks in the light profiles. The chemical abun- 
dance analysis of these outer disk regions therefore provides 
essential information on the processes that lead to the forma- 
tion of galactic disks. 

We have studied the oxygen abundance gradient in the outer 
disk of M83, in the radial span between 0.6 and 2.6 times the 
isophotal radius, obtaining deep optical spectra for nearly 50 
Hll regions. We have obtained nebular O/H ratios adopting 
different strong line metallicity indicators, as well as the direct 
method that involved the use of the [O III] A4363 auroral line. 
Regardless of the abundance indicator used, we obtain a flat 
oxygen abundance gradient beyond the R25 isophotal radius, 
and a drop in abundance slightly beyond this galactocentric 
distance. The value of the approximately constant abundance 
in the outer disk depends on the adopted diagnostic, and varies 
between 12 + log(0/H) = 8.2 and 12 + log(0/H) = 8.6. 

By analogy with the results obtained for LSBGs, we can 
interpret the observed abundances in the outer disk of M83 as 
due to its relatively unevolved state. The slow evolution of 
the outer disk, perhaps still in the process of being assembled 
via gas inflow, results in a flat distribution of chemical abun- 
dances outside the region of main star formation. The rather 
abrupt change in star formation properties and gas density be- 
tween the inner and outer disk produces the observed break 
near the optical edge of the galaxy. 

Generalizing the results obtained for M83 is, at this point, 
obviously premature. A conclusion on whether the trends ob- 
served in M83 are a feature of the majority of extended spiral 
disks, or are instead peculiar to this galaxy, awaits further de- 
tailed chemical abundance studies in a larger number of simi- 
lar systems. 
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APPENDIX 

We summarize in Table|4]the equatorial positions (Columns 2 and 3) and the redshifts (Column 4) of background emission-line 
galaxies that were serendipitously discovered in this study. The numbering scheme for the FORS field in Column 5 is the same 
as in Fig.Q] 
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